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Abstract 

The Constant Elasticity of Variance (CEV) model is mathematically 
presented and then used in a Credit-Equity hybrid framework. Next, we 
^ ' propose extensions to the CEV model with default: firstly by adding a 

l/^ , stochastic volatility diffusion uncorrelated from the stock price process, 

^D ' then by more generally time changing Bessel processes and finally by 

C^ I correlating stochastic volatility moves to the stock ones. Properties about 

"N] , strict local and true martingales in this study are discussed. Analytical 

7"^ ' formulas are provided and Fourier and Laplace transform techniques can 

j5^ I then be used to compute option prices and probabilities of default. 

^ ; 1 Introduction 

s 



X 



It has been widely recognized for at least a decade that the option pricing 
^ . theory of Black and Scholes (1973) and Merton (1973) is not consistent with 

market option prices and underlying dynamics. It has been noted that options 

with different strikes and maturities have different implied volatilities. Indeed, 

C^ ' markets take into account in option prices the presence of skewness and kurtosis 

in the probability distributions of log returns. In order to deal with those 
effects, one could use stochastic volatility models (e.g. Heston (1993), Hull 
and White (1987) or Scott (1987)). Another common alternative is to use a 
deterministic time and stock price dependent volatility function, the so-called 
local volatility to capture these effects. One would then build the volatility 
surface by excerpting the values of this function from option prices, thanks to 
the well-known Derman and Kani (1994) and Dupire (1994) formula. 

*We thank Marc Yor for all the important contributions to this paper. Any remaining 
errors are our own. 



One of the first models developed after Black Merton Scholes (1973) is the 
Constant Elasticity of Variance model pioneered by Cox (1975) where the volatil- 
ity is a deterministic function of the spot level; This latter model is somehow 
an ancestor of local volatility models. It has very interesting features since it 
suggests that common stock returns are heteroscedastic and that volatilities im- 
plied by the Black and Scholes formula are not constant, in other words skew 
exists in this model. Another interesting property is that it takes into account 
the so called "Leverage Effect" which considers the effects of financial leverage 
on the variance of a stock: a stock price increase reduces the debt-equity ratio 
of a firm and therefore decreases the variance of the stock's returns (see for 
instance Black (1976), Christie (1982) or Schwert (1989)). A last but not least 
feature of this model is that it has a non-zero probability of hitting and this 
could be of importance when one is interested in modelling default by defining 
bankruptcy as the stock price faUing to 0. 

For the last few years, the credit derivatives market has become more and 
more important and the issue of modeling default has grown, giving birth to 
two main classes of models. The first class is the structural models of the firm 
pioneered by Merton (1974) where bankruptcy occurs if the asset value falls to 
a boundary determined by outstanding Habilities. Other early work on such 
models was done by Black and Cox (1976) and Geske (1977). The other class 
commonly called reduced-form models is less ambitious than structural models. 
They consider the time of default as an exogenous parameter that they calibrate 
under a risk neutral probability to market data. These models were developed 
by Artzner and Delbaen (1995), Jarrow and Turnbull (1995), Duffie, Schroder 
and Skiadas (1996) and Madan and Unal (1998). 

The credit risk is also a component of the equity derivatives market as it may 
appear in convertible bonds or more generally in Capital Structure Arbitrage 
for people that embedded it from out-of-the money puts. It is then clear that 
having a consistent modeling of equity and credit is essential to eventually be 
able to manage those cross-asset positions. Indeed, a market standard has been 
developed during the last few years which involves a jump diffusion dynamics 
for the stock price with a local probability of default for the jump factor. This 
kind of model has been presented for instance in Ayache, Forsyth and Vetzal 
(2003). An important drawback of this modeling is that the stock has to jump 
to zero in order to default, which isn't a realistic assumption as we can see on 
several historical data and as argued in Atlan and Leblanc (2005). 

The necessity to have stock price diffusions that don't jump to zero in order 
to default and still have a non-zero probability of falling to zero leads us to nat- 
urally consider CEV processes. Moreover, CEV models have the advantage to 
provide closed-form formulas for European vanilla options and for the probabil- 
ity of default. Those computations were originally performed by Cox (1975) in 
the case where the stock can default and by Emanuel and McBeth (1982) when 
the stock never defaults. Then, one may want to add a stochastic volatility 
process to the CEV diffusion in order to capture some volatility features such 
as a smile or such as a more realistic volatility term structure. Finally, to get 
more dependency between the stock price and the volatility, one may add some 



correlation. 

Those guidelines lead us to study in section 2 the one-dimensional marginals, 
the first-passage times below boundaries and the default of martingality of Con- 
stant Elasticity of Variance processes, mainly by relating those latest to Bessel 
processes. In section 3, we propose a CEV model that is stopped at its de- 
fault time and we provide closed form formulas for European vanilla options. 
Credit Default Swaps and Equity Default Swaps. Section 4 extends the Con- 
stant Elasticity of Variance framework to a Constant Elasticity of Stochastic 
Variance one by firstly adding a stochastic volatility to the CEV diffusion and 
in a second time more generally consider time-changed Bessel processes with a 
stochastic integrated time change. Quasi-analytical formulas conditionally on 
the knowledge of the law of the time change are provided for vanilla options and 
CDSs and examples are given. Section 5 adds a correlation term to the general 
time-changed power of Bessel process framework, once again quasi analytical 
formulas conditionally on the knowledge of the joint law of the time change 
and a process related to the rate of time change are provided for probabilities 
of default and for vanilla options, and computations for several examples are 
shown. All the models proposed in this paper are true martingales and the 
martingale property is carefully proven for the different frameworks. Finally, 
section 6 concludes and presents possible extensions of this work. 

Convention For strictly negative dimensions we define squared Bessel processes 
up to their first hitting time of after which they remain at 0. 

We set this convention because we wish to consider positive Bessel processes. 
For a study of negative dimension Bessel processes with negative values, we refer 
to Going- Jaeschke and Yor (2003). 

2 A Mathematical Study of CEV Processes 

2.1 Space and Time Transformations 

A reason why Bessel processes play a large role in financial mathematics is that 
they are closely related to widely used models such as Cox, Ingersoll and Ross 
(1985), i.e. the CIR family of diffusions for interest rates framework, such as the 
Heston (1993) stochastic volatility model or even to the Constant Elasticity of 
Variance model of Cox (1976). They are more generally related to exponential 
of time-changed Brownian motions thanks to Lamperti (1972) representations. 

Let us now concentrate on the CIR family of diffusions: they solve the 
following type of stochastic differential equations: 



dXt = {a- bXt)dt + a^\x7\dWt (1) 

with Xo = xq > 0, a £ IR+, & e M, cr > and Wt a standard Brownian motion. 
This equation admits a strong (e.g. adapted to the natural filtration of Wt) 



unique solution that takes values in R+. 

Let us remark that squared Bessel processes of dimension S > can be seen 
as a particular case of a CIR process with a — 6, b = and cr = 2. We also 
recall that a Bessel process Rt solves the following diffusion equation 

A 1 

dRt = dWt + -T^dt 
I tit 

where for 5 — \, the latter ^^dt must be replaced by a local time term. 

One is now interested in the representation of a CIR process in terms of a 
time-space transformation of a Bessel Process: 

Lemma 2.1 A CIR Process Xt which solves equation 0) can he represented in 
the following form: 

Xt^e-''BESQ^s,..){^{e''-l)) (2) 

where BESQ(^s,xo) denotes a squared Bessel Process starting from xq at time 
t — of dimension (5 = || 

Proof : This lemma results from the identification of two continuous functions 
/ and g (with g strictly increasing and g{0) = 0) such as 

Xt^fit)BESQ^s.,„~,{g{t)) 

To do so, we apply Ito's formula and Dambis (1965), Dubins-Schwarz (1965) 
theorem ■ 

This relation is widely used in finance, for instance in Geman and Yor (1993) 
or Delbaen and Shirakawa (2002). 

Let us now introduce the commonly called CEV (Constant Elasticity of Vari- 
ance), which was introduced by Cox (1975, 1996) and that solves the following 
equation: 

dXt = ^iXtdt + aX^dWt (3) 

with Xq = a;o > 0, a e M, /i G M, ct > and Wt a standard brownian motion. 

Lemma 2.2 A CEV Process Xt which solves equation \^ can he represented 
as a power of a CIR process, indeed for [3 = 2{a — 1), 1/ X^ solves 



d{^) = {a-b^)dt + j:J\^\dWt (4) 

where a = '^ ^ ''' , b — /3/i, S = — /3cr and . 

Proof : This lemma is just an application of Ito's Lemma. ■ 

As a consequence of Lemma 2.1 and Lemma 2.2, one obtains the following 
representation for a CEV process: 



Proposition 2.3 A CEV Process Xt solution of equation ^ can be represented 
in the following form: 

where BESQ^s,xo) denotes a squared Bessel Process starting from xq at time 
t = of dimension S. 

2.2 Distributions and Boundaries 

We will now recall well known results about squared Bessel processes and deduce 
some properties about CEV processes. 

Path Properties 

Proposition 2.4 According to its dimension, the squared Bessel process has 

different properties: 

(i) if S < 0, is an absorbing point. 

(a) if S < 2, {0} is reached a.s. 

(Hi) if S > 2, {0} is polar. 

(iv) if d <2, BESQ is recurrent. 

(v) if S > 2, BESQ is transient. 

(vi) ifO<S<2, {0} is instantaneously reflecting. 

Proof : The proof can be found in Revuz and Yor (2001). ■ 

As a consequence, one may give some properties of the CEV diffusions. A 
topic of interest for the remaining of the paper is whether or not {0} is reached 
by a CEV process. 

Proposition 2.5 According to the value of a, the CEV diffusion has different 

properties: 

(i) if a <1, {0} is reached a.s. 

(a) if a < \, {0} is instantaneously reflecting. 

(Hi) if ^ < a <1, {0} is an absorbing point. 

(iv) if a> 1, {0} is an unreachable boundary. 

Proof : It is a consequence of the previous proposition and of Proposition 2.3. ■ 



Distributional Properties 

It is important to notice that the law of a squared Bessel process can be seen 
in terms of non-central chi-square density: 



Lemma 2.6 For any BESQs,x, one has: 

BESQsAt) = tT/(^^f ) (6) 

where V^°''^' is a non-central chi-square r.v. with a degrees of freedom and non- 
centrality parameter b> 0. Its density is given by: 

Proof : This proof results from simple properties of Laplace transforms and 
can be found for instance in Delbaen and Shirakawa (2002). ■ 

We leave to the reader the calculation of the CEV density in terms of non- 
central chi-square distributions. 

Let us recall a useful result for the remaining of the paper on the moments 
of a squared Bessel process: 

Corollary 2.7 // y ("■''•' is a non-central chi-square r.v. with a degrees of free- 
dom and noncentrality parameter b > 0, then for any real constants c and d: 

niy^-'^niv...>_.,] - e-h^ Y. {'^r^-}0^Gin + ^ + c, ^) (8) 

n>0 V 2 '' 

where G is defined as follows: 

Proof : This calculation is a simple application of Lemma 2.6. ■ 

Finally, for the computations involved in this paper, one recalls the two 
following identities on the complementary non-central chi-square distribution 
function Q that one can find in Johnson and Kotz (1970): 

0(2z, 2i^,2K) == Yg{n,K)G{n + i^ -l,z) 

ri>l 

1- Q(2k, 2z^-2,2z) = ^5(n + z/- l,K)G(n,z) 



n>l 



where g{x,y) = -^{x,y). 

First-Hitting Times 

We now concentrate on the first hitting time of by a Bessel process. For this 
purpose, let us consider a Bessel Process R of index j^ > starting from at 
time 0, then, one has: 

LAR) '=' ^ (9) 



where Li{R) — sup{t > 0,Rt — 1} and Z^, is a gamma variable with index v 
that has the following density: 

nz, e dt) = ^ wj) Mt>o}dt (10) 

This result is due to Getoor (1979). Thanks to results on time reversal (see 
Williams (1974), Pitman and Yor (1980) and Sharpe (1980)), we have: 

[Rto-u] u < ToiR)) '=^ (i?„; u < Li{R)) (11) 

where -R is a Bessel Process, starting from 1 at time of dimension S — 2{1 ~ ly) 
and Tq{R) = inf{t > 0, i?t = 0}. As a consequence, one has: 

UR) ^^ ^ (12) 

Using the scaling property of the Squared Bessel Process, one may write: 

MBESQi) ^^ ^ (13) 

with (5 = 2(1 - ly). 
Hence, we are now able to state the proposition below: 

Proposition 2.8 The probability of a CEV diffusion solution of equation l^ 
to reach at time T with a < I is given by: 

P(To < T\X„ = xo) = G(^^^^, Ct) (14) 

where G and £,t are defined as follows: 

G{x,y) = / \z>o}dz (15) 

Jz>y ^ \^) 

2(l-a) 

^^ (l-a)a2(l-e2("-i)MT) '^ ''> 

Proof : This proof is just a consequence of Proposition 2.3 and equation ((T3|l . 



Remark 2.9 The calculation of the probability of default was originally done 
by Cox (1975). 

In order to compute first-passage times of scalar Markovian diffusions below 
a fixed level, let us recall Ito and McKean (1974) results. If {Xt,t > 0) is scalar 



Markovian time-homogeneous diffusion starting from xq at time of infinitesi- 
mal generator C and that we define th — inf{i > 0;Xt < H} for H < xq, then 
for any A > 0, we have 

where (fix is solution of the ODE 

C(f>^X(l) 

with the following limit conditions: 

limx^ac 4>\{x) = 

If is a refiecting boundary then 4>\{0+) < oo 

If is an absorbing boundary then (/)a(0+) = oo 



As a first example, let us now consider the first-hitting time below a fixed 
level < y < X of a Bessel process Rt of dimension S = 2{v + 1) starting from 
x : 

Ty = inf{t >0:Rt< y} 

The law of Ty (see It6 and McKean (1974), Kent (1978) or Pitman and Yor 
(1980)) is obtainable from the knowledge of its Laplace transform C. One has 
for any positive A 

£{X) = E[e-^^«] 

x-^K^ixV^J^) 



y--K^iyV2X) 
where f £ M\Z and K,y is a Modified Bessel function defined as follows: 



kW{iy + k + l) 



As a second example that will be useful for the computation of EDS prices, 
let us write the infinitesimal generator of a CEV process: 

LcEv4>^(yx -r^ + tJ'X-- 
ax^ ax 

that must solve 

CcEv4> — A0 

with the following conditions: 

4>\ is a decreasing function, limj;^oo 4>\{x) = 



If a < ^, then (f)\{0+) = oo 
If ^ < a < 1, then 0a(O+) < oo 



We obtain the following result whose computations of the Laplace transforms 
were originally performed by Davydov and Linetsky (2001): 

Proposition 2.10 For a CEV process solution of |0) with a < 1 and ^ ^ 0, 
then 

MA = .-* exp ( - ^^^ »■,„ (l*f^) (17) 



where 



fc = .gn(,)(,,A^-l)-., ,/ . - 



^4(1 -a) 2' 2\ijL\{l~a) 4(1 - a) 

and Wk,m is a Whittaker function 

The definition of the Whittaker function can be found for instance in Abramowitz 
and Stegun (1972). 

2.3 Loss of Martingality 

Let us now state a result on some martingale properties of Bessel processes 
which play an essential role in pricing theory as is well known: 

Theorem 2.11 Let Rt be a Bessel process of dimension S starting from a ^ 0, 
then: 

(i) If S <Q, Ri" is a true martingale up to the first hitting time ofO. 
(a) IfO<S<2, the process i?j^ ~ Lt is a martingale where Lt is a continuous 
increasing process carried by the zeros of {Rt, t > 0). 
(Hi) If 5 — 2, log{Rt) is a strict local martingale. 

(iv) If 5 > 2, R^^ is a strict local martingale. Moreover, the default of mar- 
tingality is 

7^'^) (i) = E[i?2-*] - E[i?2-'5] = a^-^¥\^-^^ (To < t) (18) 

where P^ is the law of {Rf\t> 0). 

Proof : (i) and (ii): Since {0} is reached a.s., we need to apply Ito's formula 
in a positive neighborhood of 0. Let us consider e > 0. We have: 



(e+i??)i-i = {eWY---+{2-5) [ {e+R^r^ RsdW^+eSil-^) [ 

Jo ^ Jo 



ds 



{e + Ri)i 



Then, as e tends to zero, it is easy to see the first term of the right hand side is 
a true martingale for S < 2 and that the second term of the right hand side is 
increasing whose support is the zeros of {Rt,t > 0) when S >0. li S < 0, T^~ 



is a true martingale. 

(iii): By applying Ito formula, we obtain 



logiRt) = log(i?o) + / 
Jo 



Rs 



We then see that log(i?t) is a local martingale. We prove that it is a strict local 
martingale by first using the fact that 

then writing that 

E[log(i?0] = ^E[logiR^t)]^^E[log{BESQ2fi{t)+BESQaAt))] 
> ^E[\og{BESQ2Am 
and finally since BESQ2,o{t) = 2te where e is a standard exponential, we obtain 

E[log{Rt)] >C+^ log(t) -^t^oo +CX3 

which shows that log(i?i) is not a true martingale, 
(iv): To compute 7' •*, we will need the following result: 

Lemma 2.12 Let {Rl ,t > 0) be a Bessel process of dimension S > 2 starting 
from a ^ 0, then 

/ p(<5) ^ 2-S 
l4-5 _ ^t \ m,S 



'' a\ntn{t<To} - [ ^ J '^alUt 

where TZt is the canonical filtration of the Bessel process and Tq the first-hitting 
time of the level 0. 

Proof : This property results from a double application of Girsanov Theorem 
by computing 

jm4-(5 ipS 

"■'^a\TZtn{t<To} , '^'^alTZt 



dp2|„ dp2,„ 

a\TCt a\TCt 

Then, by identification, one gets the announced result. A more general result 
can be found in Yor (1992). ■ 

We may then write 

and consequently compute the default of martingality. ■ 

A proof in the case < 6 < 2 can be found in Donati-Martin et al. (2006) 
and proofs when S > 2 exist in Elworthy, Li and Yor (1999). As a consequence, 
we obtain similar results for a CEV process. 



10 



Proposition 2.13 Let Xt be a CEV Process of elasticity a solving the following 
equation 

dXt = tidt + aX^dWt 

then: 

(i) If a < \, the process e~^*Xt is a true martingale up to the first hitting time 

ofO. 

(a) If -^ < a < 1, the process e ^^Xt — L-^ is a martingale where Lf is a 

continuous increasing process carried by the zeros of {Xt,t > 0) and consequently 

e"^*X( is a true martingale up to the first hitting time ofO. 

(a) If a = I, e^^^Xt is a geometric Brownian motion and hence a martingale. 

(Hi) If a > 1, e^''*Xf is a strictly local martingale. Moreover, the default of 

martingality is 

^x{t) - E[Xo] - Eie-f^'Xt] = a;oG(— ^— ,Ct) (19) 

Z[a ~ ij 

where G and Ct are defined as follows: 



G{x,y) = / l{2>o}C^^ 

Jz>y -L \^) 

2(1-q) 

Ct ~ 



(a-l)a2(e2("-i)MT_i) 

Proof : This is just an application of Theorem 2.11, Proposition 2.3 and equa- 
tion Jul. ■ 

A proof of the failure of the martingale property can be found in Lewis (1998). 

Remark 2.14 For a > I, one has V(t, K) e M^; 

E[(e"'^*Xt - K)+\ - n{K - e-^*Xt)+] + lx{t) = E[Xo] - K (20) 

The last equation shows that in the case of a strictly local martingale, the 
Call price must incorporate the default of martingality in order to remain in a 
No Arbitrage model. For a study on option pricing for strict local martingales, 
we refer to Madan and Yor (2006) for continuous processes and to Chybiryakov 
(2006) for jump-diffusion processes. Lewis (2000) also did this study in the 
case of explosions with stochastic volatility models and in particular for a CEV 
diffusion. 

3 Credit-Equity Modelling 

3.1 Model Implementation 

Usually, in the mathematical finance literature, one defines a CEV diffusion for 
the stock price dynamics S to be 

—i = ^dt + crS^-^dWt 
11 



First of all, in a credit perspective, we will just consider the case a < 1 since we 
are interested in models with a non-zero probability of default. Once the stock 
has reached zero, the firm has bankrupted and that is the reason why we stop 
the CEV diffusion at its first default time. Then from what has been proven 
above, we know that the stock price process hence defined is a true martingale 
and that ensures the Absence of Arbitrage and moreover the uniqueness of the 
solution. Hence, the stock price diffusion now becomes under the risk-neutral 
pricing measure: 

-J.^rdt + aS^-^dWt if t < t. 

St^O if t>T. 

where r = To (5) — inf{t > 0,S't = 0}. In other words, the stock price process 
considered is nothing else than a stopped CEV diffusion {StAT)t>o- 

Remark 3.1 Delbaen and Shirakawa (2002) showed the existence of a risk- 
neutral probability measure whose uniqueness is only ensured on the stock price 
filtration considered at time t Tr = a{St,t < r). Since our purpose is to 
compute the price of options whose payoffs are Tt — measurable, we have the 
uniqueness of the no-arbitrage probability. 

3.2 European Vanilla Option Pricing 

Lemma 2.12 states that 

K\ntn{t<n} = \^) ' '^^l^' ^^^^ 

Thanks to this identity, we obtain the law of the stopped CEV diffusion at a 
given time. Lemma 2.6 and Corollary 2.7 enable us to compute the call and put 
option price: 

For the call Co option price 



Co = e-^'^n{STAr-K)+\ 
= e-^^niSr - K)+lT<r\ 



and the put Po option price: 

PO = e-^-'^niK - STAr)+\ 

= e-"^E[(i^ - StJ+Itkt] + Ke-'-^FiT < T) 
Consequently, for the call price: 

Co - SoQ{zT, 2 + -^, 2Ct) - Ke-'-^il - Q(2Ct, -^, zt)) 
i — a 1 — a 



12 



and for the 


put price: 












R, 


= Ke- 


-^^(Q(2Ct, 


1 


-,zt) 
a 


^^2{l-a) 


Xt)) 




J^o 


1- 






-So 


(1 - Q{ZT, 


2 + 


1 


2Ct)) + Ke-^' 


^P(r < T) 




1-a' 




= Ke- 


■^^0(2Ct,- 


1 


-iZt) - 


-So(l-Q(zT. 


.2 + - 


1 



a 



,2Ct)) 



where 



2ri^2(i- 



ZT 



Ct 



a2(l-a)(e2(i-")'-^-l) 

„c2(l-a) 

^jfo 

(l-a)(T2(l_e-2(l-a)rT) 



Hence, one easily verifies that the put-call parity is satisfied. Closed-form 
CEV option pricing formulas were originally computed by Cox (1975) for a < 1 
and Schroder (1989) expressed those formulas in terms of non-central chi-square 
distributions. Computing option prices using the squared Bessel processes dis- 
tributions was done by Delbaen and Shirakawa (2002). 

3.3 Pricing of Credit and Equity Default Swaps 

Since we are dealing with default probabilities, it is obvious to consider deriva- 
tive products relying on these probabilities. One of the most liquid protection 
instruments against default is the Credit Default Swap (CDS). The buyer of the 
protection agrees to pay periodical amounts until a default time (if it occurs) 
and in exchange receives a cash amount which is a notional amount minus a re- 
covery rate in the case the company on which the contract is written, defaults. 
The payoff of such kind of contract is: 



Hcds = - V e-'-^'Cl{,>T,} + e---il - i?)l{,<T„} 






where C is the periodical coupon, Ti, ...,T„ the payment dates, R the recovery 
rate assumed to be deterministic and r the default time. For simplicity purposes, 
we consider in this paper deterministic interest rates. The CDS Fair Price is 
the expectation of the payoff conditionally to the spot price filtration taken at 
the pricing time, e.g.: 

n 

By absence of arbitrage, one must have CDSt{Ti,Tn; C; R) — and then 

(l-i?)E[e--(--*)l{,<T„}|5,] 



C 
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From Proposition 2.8, we know the value of (P(t > Ti\St)i<i<n)- It then remains 
to compute the following quantity E[e~'"'^li-<t] to be able to price the CDS 
coupon C. By an integration by parts, we show that 



ft 

-rtji 



E[e"'^^l^<t] = e-''*P(T <t)+r e^'^"P(T < s)ds (22) 

Otherwise, one could just obtain this expectation by directly using the density of 
the first-hitting time of that is provided by the differentiation of the cumulative 
distribution function : 



^ 2r(l-a)C/'^-°'e-C^ 

where Q is defined above. 

EDSs are very similar to CDSs except that payouts occur when the stock price 
falls under a pre-defined level, which is often referred to as a trigger price. The 
trigger price is generally between 30% and 50% of the equity stock price at the 
beginning of the contract. Hence, these contracts provide a protection against a 
credit event happening on the equity market for the buyer. They were initiated 
by the end of 2003. At that time, it had become difficult in many countries to 
structure investment-grade credit portfolios with good returns because the CDS 
spreads were tightening, as reported by Sawyer (2003). Another reason why 
people have interest in those contracts is because the settlement of the default 
is directly observed on the stock price. Let us now define tl as the first passage 
time of the stock price process under the level L < Sq. Formally, we write 
Tl — inf{t > 0; S't < L}. We recall the general valuation formula of an EDS: 

n 

EDStiT,,T^;C;R) - -C^ e-'^(^'-*)p(Ti > T.l^,) + E[e-'^^"""*'l{r,<T„}|5t] 

where C is the coupon, Ti, ...,T„ the payment dates and r the risk-free interest 
rate. Again, by absence of arbitrage, we can find the coupon price, by stating 
that at the initiation of the contract: 



Or equivalently 



C 



EDSt=o{Ti,Tn;C;R)^0 



In order to price the coupon C, one needs to evaluate: 

E[e-''^^l{,^<t}] and P(tl < t) 
An integration by parts gives the Laplace transform of P(tl < t) for any A > 



/ dte-^^PiTL <t) = 
Jo 
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Applying Fubini theorem, one observes that 



+00 Efp"(''+^)'^il 



HrL<t}\ 



A 

Hence using Proposition 2.10, one is able to compute the Laplace transform 
of the desired quantities necessary to evaluate an EDS. One can then use numer- 
ical techniques (see Abate and Whitt (1995) for instance) to inverse the Laplace 
transform in order to evaluate prices. 

4 Stochastic Volatility for CEV Processes 

4.1 A Zero Correlation Pricing Framework 
Impact of a Stochastic Time Change 

Due to the very important dependency between the probability of default, the 
level of volatility and the skewness, we were naturally brought to consider ex- 
tensions of the CEV model that could relax the high correlation between these 
three effects. More precisely, in a CEV model, if one first calibrates the implied 
at-the-money volatility, then either the skewness or the CDS will be calibrated 
on adjusting the elasticity parameter. Hence, to be able to get some freedom on 
the volatility surface, a possible extension is to introduce a stochastic volatility 
in the CEV model instead of a constant volatility. A CEV diffusion with a 
stochastic volatility is actually just a power of a squared Bessel Process with a 
stochastic time change instead of having a deterministic one like in Proposition 
2.3. 

Another extension is to consider a power of a Bessel Process time changed 
by an independent increasing process. More precisely, one writes the following 
process for the stock price: 

St^e^'BESQ'-j^i^t) if t<T. (23) 

St=0 if t> T. 

where x = sf^ , t = To{S) = M{t >0,St = 0} = ^-\To{BESQ)) and ^ is an 
strictly increasing continuous integrable process independent from the squared 
Bessel process. Subordinating a continuous process by an independent Levy 
process is an idea that goes back to Clark (1973). Stochastic time changes are 
somehow equivalent to adding a stochastic volatility in stock price diffusions. 
The basic intuition underlying this approach could be foreseen through the 
scaHng property of the Brownian motion, or through Dambis (1965),Dubins 
and Schwarz (1965) (DDS) theorem or even its extension to semimartingales 
by Monroe (1978). More recently, Carr et al. (2003) generated uncertainty 
by speeding up or slowing down the rate at which time passes with a Levy 
process. Our approach differs from the one done in the Levy processes literature 
for mathematical finance: We are not considering the exponential of a time 
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changed Levy process but a power of a time changed Bessel process. Thanks 
to Lamperti representation (1972), this means that we are considering a time 
changed geometric Brownian motion B. More precisely, it is known that 

/■* ds 
Rt = exp (Bct + vCt) and ^t = ~52 

Jo Rs 

where {Rt,t > 0) is a Bessel process of dimension S = 2(1 + u) starting from 
a ^ 0. Hence the time change considered in the stock price is 

f^' ds 



/O ^s 

and the stock price process as defined in equation i2'A\i can be identified as 
follows: 

St^e^'exp{^2iyBY,-^Yt) 

As a consequence, we have now proposed a new class of time changes where 
analytical computations are possible thanks to a good knowledge of Bessel pro- 
cesses. 

For the absence of arbitrage property, there must exist a probability under which 
all the actualized stock prices are martingales. A very simple property on mar- 
tingales is that a process Mt is a martingale if and only if for every bounded 
stopping time T, E[Mt] = E[Mo]. Nonetheless, this result is not very conve- 
nient. Let us state and give a straightforward proof of the martingality of the 
stock price process 

1-- 
Proposition 4.1 Consider Mt — BESQ ,^ \{S^t/\T) where following the previ- 
ous hypotheses £,t is a strictly increasing continuous integrable process indepen- 
dent from BESQ, T is the {Mt,t > 0) first hitting time of and BESQ(^s,x) ** 
a squared Bessel process of dimension S starting from x ^ Q, then {Mt,t > 0) is 
a true martingale. 

Proof : Let us define TZt — cr{Rs; s < t). We then naturally write the canon- 
ical filtrations TZ^^ — cr(i?{js < t) and St — a{^s', s < t). For any bounded 
functional F, we want to compute 

nF{R^^-u<s){Rl;'-Rl;')] 

Since ^ is integrable and independent from R, we obtain by using Fubini theorem 
E[F(i?5„ ■u<s) {Rj-^ - i?!;"^)] = E E[F(i?^„ ;u<s) (i?^-'' - i?^-'') | Et 



PH,(da)E[F(i?,(„);« < .) (i?^-^^ - i?f-^^)] 



The latest quantity is null by Theorem 2.11 and we have then shown that for 
s <t < T 

Rl;^ = E[Rl-^\n^^] 

which is the announced result. ■ 
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Pricing Vanilla Options 

One can find closed-form formulas for the call and put options prices. Let us 
define the two following quantities Co(x, S, K, T; So) and Po{x, S, K, T; Sq): 

Co{x,5,K,T;So) = SoQC- '- ,4 -,5,^^) 

X X 

-Ke-^^{1 - Q{^^^,2 - S, ^— > )) 

X X 

_2 2 

Poix,S,K,T;So) = Ke-'''^Q{^^,2 - d }^^ ^ ^"' ) 

X X 

-So{i-Q r' ' A-S,^) 

X X 

From there, one may obtain the option prices under the new general framework. 
Proposition 4.2 // one has the following stock price process: 

St^e'-'BESQ\ll^{it) if t<T. 
St=0 if t> T. 

2 

where x = Sq^'' , t — Tq{S) — inf{t > 0, St = 0} and ^t is a strictly increasing 
continuous integrable process independent from BESQ whose probability mea- 
sure is /ijj(dx), then: 



Co = / Co{x,S,K,T;So))i^^{dx) 
Po = / Po{x,S,K,T;So)^J.^Tidx) 



Proof : Let us prove this result for the call option price, a similar result may 
be obtained for the put price. One has: 

Co = e-'-^EiiST - K)+] 

= e-'-^E(E[(^T-if) + k(6;5<T)]) 
= E[Coi^T,S,K,T;So)] 



Computing the Default 

Having the integrability of the change of time and knowing its density, one could 
find a closed-form formula for the probability of default r — To{S) — mi{t > 

0,St = 0} where ^4 ~ e^^BESQ,g ^s{£,t). Let us now compute the probability 
of default the proof of which is left to the reader: 
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Proposition 4.3 // one considers a stock price process defined as follows: 
then the probability of default t — inf {t > 0,St — 0} is given by 

s sf" 



I 1;t 



P(t < T) = E 
where G is the complementary Gamma function. 

4.2 CESV Models 

Stochastic volatility models were used in a Black and Scholes (1973) and Merton 
(1973) framework mainly to capture skewness and kurtosis effects, or in terms 
of implied volatility skew and smile. In a Constant Elasticity of Variance frame- 
work, one would use stochastic volatility not to capture the leverage effect which 
partly already exists due to the elasticity parameter but to obtain environments 
for instance of low volatilities, high probabilities of default and low skew. Let us 
consider an integrable jump-diffusion process (ut, t > 0) to model the volatility. 
We will call those diffusions Constant Elasticity of Stochastic Variance (CESV) 
for the remainder of the paper. Leblanc (1997) introduced stochastic volatility 
for CEV processes. 

Hence, the class of models under a risk-neutral probability measure proposed 
is of the following form: 

^ = rdt + at-Sr'dWt 
St 

where a is assumed to be independent from the Brownian motion driving the 
stock price returns. Next, within an equity subject to bankruptcy framework, 
we are going to stop the diffusion when the stock reaches just as in the previous 
section. As a consequence, our diffusion becomes: 

—^ = rdt + at-S^-^dWt if t<T. 

St = if t>T. 

where r = To(S') = inf{i > 0, 5* = 0}. 

Before giving any concrete examples, let us show how CESV models can be 
seen as Bessel processes with a stochastic time change. So first, let us recall 
elementary results: 

Lemma 4.4 Let R be a time change with s y^ Rg continuous, strictly increas- 
ing, Rq = and Rt < oo, for each t > 0, then for any continuous semimartingale 
X and any caglad (left continuous with right limits) bounded adapted process H , 
one has: 

H,dX, = / HrJXr^ (24) 

Jo 
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Proof : The proof can be found in Revuz and Yor (2001). ■ 

Then, using Lemma 4.4, (DDS) theorem and Ito formula, we obtain that 



St = e^'^BESQ^''-"^ ,,1 „, (HtAr) 

^(2-l/(l-a),S^<'-°')^ '^'^^^ 

T = inf {t >0,St = 0} 

Ht = {1-af f a^^e-^^^-^>'ds 
Jo 

Ht is by construction an increasing continuous integrable process. 
Hence {e~^*St,t > 0) is a continuous martingale by Proposition 4.1. All the 
results of the previous subsection apply and we are able to compute Vanilla 
option and CDS prices conditionally on the knowledge of the law of Ht. As a 
result, we showed that a CESV model is in fact a timed-changed power of Bessel 
process where the subordinator is an integrated time change Ht ~ J^ hgds with 
a specific rate of time change ht that is defined by 

ht = {l- a)2a2e-2(i-")'-* 



We now provide two examples of well-known stochastic volatility models 
where we compute the law of the time change. 

Heston Model Let us first consider a CIR (1985) difi^usion for the volatility 
process 

da^ = k{0 -~ at)dt + -qatdW" and al ^ x > 

where k, 9 and 77 are strictly positive constants and W^ is a Brownian motion 
independent from W . In fact we are proposing a variation of the Heston (1993) 
model by considering a ^ 1. We then want to compute the law of 

Tj (A \2 I 2 -2(1-Q)rs 7 

Ht — [1 — a) / CT^e ^ ' ds 
Jo 

More precisely, we will compute its Laplace transform, that is to say, for any 
A>0 

For this purpose, let us use the following result: 

Lemma 4.5 If X a squared Bessel process BESQ(s.x) starting from a; 7^ and 
of dimension 5, then for any function f : IR+ -^ M+ such that for any t > 0: 
/o f(s)ds < 00, we ho 



we nave 



E 



exp ( - I' X.m,.)] = -Vjij^ exp I (*',(0) - ||^) 
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where (j)f is the unique solution of the Sturm- Liouville equation 

<j,"f{s) = 2f{s)cj,f{s) 
where s G [0;oo[, 0/(0) — 1, 4>f is positive and non-increasing and 

/■* ds 



Proof : The proof can be found in Pitman and Yor (1982). 
By Lemma 2.1, we can see that 

du 



V Jo 



(^ + 1) 



2[-^i^ + l] 



where X is a BESQri^o ,^\. Hence for any A > 0, 



E[e- 



-\Ht 



r-l{t 



with 



'^^.Kt 



cxp( - / Xsf\{s)ds) 
Jo 



/.(t) = A(^ii^)^(!^ + l)-[^^+^l 



;(t) = ^(e«*-l) and /,(t) = A(^l^-^)'(^ + l)"-^-- (25) 

Defining a — 8((1 — a)/ri)^, b ~ '^KJ-rf' and n = —2{ ^^'^ + 1) and using 
Lemma 4.5 we are brought to the resolution of the following ordinary differential 
equation 

(/)"(x) - aX{bx + l)"(/)(x) = 

Then under the boundary conditions, one obtains (see Polyanin and Zaitsev 
(2003)): 

5I^^i/(«+2) {^{bx + l)("+^)/2) 
(Pxixj = Vbx + 1 



-'-l/(«+2)(b(«+2)J 

;^^ ^l/(»+2)(^(fex+l)("+^)/^) 
-Vbx + 1 ^ ■;= 



2VaX 



(26) 



M^) - CiVb^TT/i/(„+2)( ,/^^, (6a; + i)("+^)/^) 



b{n + 2) 
2VaA 



+c,Vb^TT ifi/(„+2)(^^f^(&x + 1)("+')/') (27) 
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where with using the fact that Il{x)K^{x) — I^{x)K'^{x) = 1/x one has 

We finally obtain the Laplace transform of Ht for any A > 

with (5=^. 

A simpler example for the forward contract is provided in Atlan and Leblanc 

(2005). 

Hull and White Model Let us now consider the Hull and White (1987) 
volatility diffusion that is driven by the following stochastic differential equation: 

where 9 and 77 are positive constants and W^ is a Brownian motion independent 
from W. Then H may be computed and after some simplifications, we obtain: 



dse^iW!+us) (28) 



If we define A^ = /^ exp2(_Bs + i^s)ds where B is a Brownian motion, we rec- 
ognize a typical quantity used for the pricing of Asian options with analytical 
formulae. Thus, we can write 

4(l-a)cr2 







V^ Jo 


where v — 


7.(9- 


- ll- 2(1 -a)r). 



and obtain its law using Yor (1992), more precisely we have V(w, v) e 



p2 



oxn (211 _ id* _ J_^ /•+°o 2x 

f^A.{u) = °P^^* * ^"^ / dxe^^'^+'^e-^^^it) (29) 



00 



where: 

t/jriv)^ [ dj/cxp(-|-)e-''™^''(^)sinh(y)sin(^) (30) 

Jo 2w V 
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4.3 Subordinated Bessel Models 

Another way to build stochastic volatiHty models is to make time stochastic. 
Geman, Madan and Yor (2001) recognize that asset prices may be viewed as 
Brownian motions subordinated by a random clock. The random clock may 
be regarded as a cumulative measure of the economic activity as said in Clark 
(1973) and as estimated in Ane and Geman (2000). The time must be an 
increasing process, thus it could either be a Levy subordinator or a time integral 
of a positive process. In this paper, we only consider the case of a time integral 
because we need the continuity of the time change in order to compute the 
first-passage time at to be able to provide analytical formulas for CDS prices. 
More generally, for the purpose of pricing path-dependent options, one needs 
the continuity of the time change in order to simulate increments of the time 
changed Bessel process. Consequently, we study the case of a time change Yt 
such as 

^t = / Vsds 



where the rate of time change {yt,t> 0) is a positive stochastic process. 

As we have seen in the previous subsection, considering a stochastic volatility 

{(Jt,t > 0) in the CEV diffusion is equivalent to the following rate of time change 



2rt 

ate 

yt 



2cS-2 



(2-5)2 

where S is the dimension of the squared Bessel process. Hence, in order to 
provide frameworks where one is able to compute the law of the time change, 
we are going to go directly through different modellings of the rate of time 
change yt- 

Integrated CIR Time change As a first example, let us consider the case 
where yt solves the following diffusion 

dyt = K{e - yt)dt + rj^dW^ 

where W^ is independent from the driving Bessel process. The Laplace trans- 
form of Yt is then defined for any A > by : 

_^y _ ^-£|i exp(-2Ayo/(K + 7Coth(7i/2))) 



Ee"^-'' = e V 



(cosh(7t/2) + ^sinh(7f/2))'^'/''' 
7 = V«^ + 27?2A 

Integrated Ornstein-Uhlenbeck Time Change We now assume the rate 
of time change to be the solution of the following SDE 

dyt = -Xytdt + dzt 
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where {zt;t > 0) is a Levy subordinator. Let ipz denote the log characteristic 
function of the subordinator Zt, then 

E[e-^1 = exp {^ay„^^^) exp (^"^ ^^^^) (^l) 

Then we can compute the characteristic function of Yt for different subor- 
dinators and we present here three examples that one can find in Carr et al. 
(2003) for which we recall below the characteristic functions: 
a) For a process with Poisson arrival rate v of positive jumps exponentially 
distributed with mean fi, we have a Levy density that is 

kz{x) = -e"~l{a;>o} 
M 

and a log characteristic function 



V'^(a;) 



1 — ZX/i 

then we obtain 

f M^dx = \og({x + -)^^^ia-Xx)^^'^^^ (32) 

b) Let us consider the first time a Brownian motion with drift i^ reaches 1. It 
is well known that this passage time follows the so-called Inverse Gaussian law 
which Levy density and log characteristic function are respectively 



V2ttx^ 



ipz{x) = f — y v^ — 2ix 



and we then get 



ibJx) , 2y/iy^ - 2ix 2\/i/2A - 2ia ( Xiv"^ -2ix) 

^ -dx = — — ^ h ..;„ arctanh'-' ^ ' 



-\x A A3/2 VV v'^\-2ic 

i^log(a — Aa;) 
A 

c) Finally, recall the Stationary Inverse Gaussian case which Levy density and 
log characteristic function are 

2 

(l + t/2a;)e- — 
^^(-) = 2V2^ '^^>°> 

I ( \ *" 

Wz\X) 



\/v^ — 2ix 
From these definitions, we obtain 



ipzix) , Vi^^ - 2ix 2ia f I Xiv"^ - 2ix) 

dx = ; ; — = arctanh ■ ' 



a-\x A A3/V/^2A-2m VV v'^\-2ia 
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5 Correlation Adjustment 

5.1 Introducing some Correlation 

We propose a time-changed Bessel process as in the previous section with some 
leverage in order to get more independence between skewness and credit spreads, 
with respect to which we add a term that contains a negative correlation (equal 
to p) component between the stock return and the volatility. Hence, let us 
consider Zt a a{hs, s <t) adapted positive integrable process such as 



¥\eP^t 



is a martingale and a general integrated time change Ht — Jq hgds such as 
E(_ffi) < oo then, we can define the stock price process as follows 



cf ^^i D r C/^2 — (5 



T = inf{t >Q;St^ 0} 
where BESQ is a squared Bessel process of dimension 6 < 2 starting from 

cl/(2-<5) 

Let us first show that the process {e~'^*'St]t > 0) hence defined is a mar- 
tingale. We know from Proposition 4.1 that BESQ^^ is a martingale. Now 
because of the independence of the processes z and BESQ 

which ensures that {e^^^S^t > 0) is a local martingale. Let us show that it is 
actually a true martingale. For this purpose, let us recall some results: 

Definition 5.1 A real valued process X is of class DL if for every a > 0, the 
family of random variables XTl{T<a} is uniformly integrable for all stopping 
times. 

We now state the following property: 

Proposition 5.2 Let Mt be a local martingale such that E|Mo| < oo and such 
that its negative part belongs to class DL. Then its negative part is a super- 
martingale. Mt is a martingale if and only j/E[Mt] = E[Mo] for all t > 0. 

Proof : The proof may be found in Elworthy, Li and Yor (1999). ■ 

All the financial assets being positive, one may use a simpler result than the 
previous property the proof of which is left to the reader: 
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Corollary 5.3 Let Mt be a positive local martingale such that IE|Mo| < oo. 
Then Mt is a supermartingale and it is a martingale if and only if E[Mi] = 
E[Mo] for allt>0. 

Consequently to prove that the actualized stock price process is a martingale 
with regards to the filtration Tt ~ Ti-Ht ^ o-{hs',s < t), we just need to show that 
for any i > 

Eie-'-'St] - So 
which is the case since 



Efe-''*S't] = E\BESQ^r7^ — r] = E E\BESQ^r7^ — .\cr(hs;s < t)] 



E[eP^ 



E[BESQJ^jJa{hs; s < t)] 






5.2 Pricing Credit and Equity Derivatives 

The computation of the probability of default is immediate from Proposition 

4.3 because 

T = inf{i > 0; 5t = 0} = inf{t > 0; BESQh, = 0} 
and then for any T > 



P(t <T)=E 
where G is the complementary Gamma function 



^ 2' Ht 



Let us compute the European vanilla option prices. For this purpose, we 
define ci'^x, y, S, K, T; So) and pI^'^x, y, S, K, T; So): 



ci'\x,y,S,K,T;So) = So 



E[eP^^] 



Q{ 



{Ke-^'''^+Py^E[eP''^])^ S^ 



4-<5, 



-Ke-'-^il - Q{^^^,2 - S, ^— ^li^^^ )) 

X X 

P^''Hx,y,S,K,T;So) = Ke-^^Q{2l^,2 - s}-^ N^J^) 

XX 



-So, 



r(i-Q( 



{Ke 



- (rT+py) jg TgP^T 1 ~) 2^r7 



A -5. 



br. 



E[e''^r] - X X 

Then, the knowledge of the joint law ^^Ht.zt for any t > enables us to compute 
the option prices as in the previous section: 



-)) 



Co 



Co (2;, y,5,K, T; So)^Ht .zt {dx, dy) 
Pq''^ {x, y, S, K, T; So)^J-Ht,zt {dx, dy) 
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5.3 Examples 

Let us go through most of the time changes presented previously and see how 
we can obtain the joint law of the couple {Ht, Zt). 

Integrated CIR Time change Let us consider the following dynamics 

dht = nid - ht)dt + ■q^/h'tdW^ 

where W^ and BESQ are independent and the stability condition ^^ > 1 is 
satisfied. Let us take 

Zt = ht + [K —)Ht 

or equivalently 

pzt = p{ho + net) +p?] ^/hldW^ " ^ / ^«^* 
Jo ^ Jo 

Hence, it is obvious that 

ppzt 



is a local martingale and it is known that it is a martingale as one may check 
using the Laplace transform below, that 

E[cxp(p77 / ^/hldW^ -^ [ f^sds)] = 1 
Jo ^ Jo 

In order to compute credit and equity derivatives prices, we then compute for 
any positive A, p the Laplace transform of 

It is well known (see Karatzas and Shreve (1991) or Lamberton and Lapeyre 
(1995)) that 



E[e-^«-M'^*] = ^, J„l ^^ ^^,^^. exp(-feoi?(t,A,A^)) 



B(t,\p) 



( cosh(7i/2) + ^^^^ sinh(7t/2))^ 
Ai(7Cosh(f ) - Ksinh(^)) + 2Asinh(^) 
7 cosh(f ) + («; + \^^) sinh(f ) 



7 = yj K^ + 2jf-\ 

Heston CESV with correlation In the same class of models, let us now 
construct z in terms of the solution of the following stochastic differential equa- 
tion 

dal ^ K{d - <jf)dt + rjCTtdWf and cr^ = x. 
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First, h is defined by 



^2g2rt/(5-2) 



(2 ~ S)^ 

Then, following the same method as in the integrated CIR time change case, 
we choose z as: 

zt = ht + [n~ j—^ - —)Ht 

Consequently, ^^p^t i is a martingale. 

Hence, it remains to evaluate for any positive t the Laplace transform of 
{Ht, Zt), that is to say for any positive A, ^ 

In order to compute the above quantity, we use the following result which ex- 
tends Lemma 4.6 that one can find in Pitman and Yor (1982). 

Lemma 5.4 If X a squared Bessel process BESQ^g r^-^ starting from x ^0 and 
of dimension S, then for any functions f and g : M+ -^ K+ such that for any 
t > 0; /p f{s)ds < CO, we have 



E 



exp(- / XJis)ds-git)Xt) 
Jo 



(^}(i) + 2g(t)V./(i))^/2 



^^^w^^ </'/W + 2.9(i)0/(i)\ 
exp-U^(O) -^ > 



2V"^'' ^'f{t) + 2g{t)^Pf{t)J 
where 4>f is the unique solution of the Sturm- Liouville equation 

cj>"f{s) = 2f{s)cj,f{s) 
where s e [0;oo[, 0/(0) — 1, 4>f is positive and non-increasing and 



i^f{t) = ^f{t) 



ds 



Taking 



(2-^)2 V4Ki + ?72 
in the above Lemma, we obtain that for any positive A, /x 



(VU^W) + 2(5^e-(«+^)VA(/(t)))*/^ 

^^ X / ^',m) + 2(^e-(-+^)V,(/(t)) X 

where noting a ~ 1^5) the functions (j)\, V'a and / are defined respectively in 
j23l,CEl and (EH. 
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Integrated Ornstein-Uhlenbeck Time Change We consider the stochas- 
tic time change Ht — /q hsds and assume that {ht',t > 0) is given by 

dht ~ —Xhtdt + dzf 

where {zt;t > 0) is a Levy subordinator. Carr et al. (2003) compute the 
characteristic function <I>(t, a, h) of {Ht,zt) for any t > and it is given by 

for any (a, b) G R^ where ipz is the log characteristic function of the subordina- 
tor. Let us first notice that 

EK^']=cxp(#,(-ip)) 

We quickly recall the computations of ^{t, a, b) for different subordinators: 
a) For a process with Poisson arrival rate v of positive jumps exponentially 
distributed with mean /x, we obtain 



i^z{x) 



-dx ^ log {{X+-) ^-'^ra + At.) U^ _|_ ;y^\ _ y^^\ A(a + At.)M + iA ] 



a + \b — \x 
b) For an Inverse Gaussian subordinator of parameter v, we have 



ipzix) 2\J v^ — 2ix 

-dx = 



a + \b — Xx A 



V^'^HM^ + Marctanh^ / ^(-^'^-) 



A3/2 Vy '^^^-^ila + Afc) 

z^log((a + \b) — Xx) 
A 

c) For the Stationary Inverse Gaussian of parameter i^, we write 

ipz (x) \/v'^ — 2ix 

-dx — 



Xb— Xx X 



^^(« + ^^) =arctanhr/ '^"^ ' ''"^ 



X^/^yiy^X-2i{a + Xb) \\] ly^ - 2i{a + Xb) 
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6 Conclusion 

Twelve continuous stochastic stock price models were built in this paper for 
equity-credit modelling purposes, all derived from the Constant Elasticity of 
Variance model, and as a consequence from Bessel processes. They all exploit 
the ability of Bessel processes to be positive, for those of dimension lower than 2 
to reach and for a certain power of a given Bessel process to be a martingale. 
We first propose to add a stochastic volatility diffusion to the CEV model, then 
more generally to stochastically time change a Bessel process in order to obtain 
a stochastic volatility effect, motivated by known arguments that go back to 
Clark (1973). Next, in order to add some correlation between the stock price 
process and the stochastic volatility, we extend our framework by multiplying 
the Bessel process by exponentials of the volatility and correcting it by its mean 
in accordance with arbitrage considerations to obtain martingale models that 
are martingales with respect to the joint filtration of the time-changed Bessel 
process and the stochastic time change itself. Hence, among the different models 
proposed based on the CEV with default model, there were first the Constant 
Elasticity of Stochastic Variance ones (CESV) taking a Hull and White (1987) 
stochastic volatility as well as a Heston (1993) one. We then proposed inte- 
grated time change models, by considering an integrated CIR time change and 
an Integrated Ornstein-Uhlenbeck time change (see Carr et al. (2003)) with 
different subordinators for the process driving the diffusion. We finally added 
correlation between stock price returns and volatilities to the models presented 
previously and provided quasi-analytical formulas for option and CDS prices for 
all of them. Let us note that we discussed the true and local strict martingale 
properties of CEV processes, that we naturally extended to the time change 
framework. 

The models presented and discussed in this paper are not specifically designed 
to cope just with Equity-Credit frameworks but they also can be used for in- 
stance for FX-rates hybrid modelling by specifying stochastic interest rates. We 
can also note that a Poisson jump to default process can be added to the CEV- 
Hke framework in order to deal with credit spreads for short-term maturities. 
Campi, Polbennikov and Sbuelz (2005) and Carr and Linetsky (2005) precisely 
considered a CEV model with deterministic volatilities and hazard rates. The 
latest paper can easily be generalized to fit in our time-changed Bessel frame- 
works. Since our goal was to concentrate on continuous diffusions, we leave the 
addition of a jump to default for further research. 
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